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Pulse formation in passively mode locked lasers is often accompanied with dispersive waves that 
form of spectral sidebands due to spatial inhomogoneities in the laser cavity. Here we present an 
explicit calculation of the amplitude, frequency, and precise shape of the sidebands accompanying a 
soliton-like pulse. We then extend the study to the global steady state of mode locked laser with a 
variable number of pulses, and present experimental results in a mode locked fiber laser that confirm 
the theory. The strong correlation between the temporal width of the sidebands and the measured 
spacing between the pulses in multipulse operation suggests that the sidebands have an important 
role in the inter-pulse interaction. 
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a. Introduction Soliton-like pulses propagating in 
nonlinear media with periodic spatial inhomogeneity can 
be a source of dispersive waves that are concentrated at 
resonant frequencies, forming spectral sidebands, some- 
times called Kelly sidebands. Several papers [ij-ll] 
studied the formation and of dispersive waves by solitons 
and their propagation in optical fibers with periodically 
spaced amplifiers. The same mechanism leads to the for- 
mation of sidebands in mode locked lasers, where inho- 
mogencities in the cavity act periodically on the pulse. 

When the power of a single pulse saturates the ab- 
sorber, the steady state of a passively mode locked laser 
tends to bifurcate into configurations where two or more 
pulses run in the laser cavity simultaneously. Since the 
early experiments that demonstrated multipulse mode 
locking 0-0] it has been observed that the pulses display 
a very rich dynamics, often forming bunches, as a conse- 
quence of complex inter-pulse interactions. The interest 
in the dispersive waves in mode locked lasers, beyond 
their prominent effect on the pulse shape, arises because 
they have often been suggested as a means of inter-pulse 
interaction in multipulse mode locked lasers @-[H| ■ Here 
we focus on the role of the sidebands in multipulse mode 
locking. 

The absorber saturation leading to multipulse steady 
states is often modeled by adding a quintic term to the 
equations of motion. In former papers fl2l [l3j we applied 
the statistical light-mode dynamics (SLD) theory to this 
system, and showed that multipulse mode locking is in 
effect a series of first order phase transitions. SLD uses 
the methods of statistical physics to analyze the dynam- 
ics of the interacting many body light mode system at an 
effective finite temperature generated by cavity noise |14l- 
ITfjj . Here we apply the SLD gain balance method [l7l| to 
derive the multipulse steady states with dispersive waves 
of mode locked lasers with cavity inhomogeneities. 

Our theoretical analysis is based on the master equa- 
tion of mode- locked soliton lasers [H, , with an ad- 
ditive noise term [2fj| . where the inhomogeneities in the 



cyclic light propagation in the cavity are modeled by a 
periodic modulation of the gain and the saturable ab- 
sorption. We first study the sidebands in a single-pulse 
steady state and show that, unlike free fiber sidebands, 
the mode locked laser sidebands reach a steady state with 
a well-defined bandwidth and a Lorentzian shape; in real 
time the dispersive waves form a wide pedestal with ex- 
ponentially decaying tails. In particular we demonstrate 
how the overall phase of each sideband depends on the 
relative phase of the gain and loss modulations. The the- 
ory is firmly supported by experimental observations in 
mode locked fiber lasers. 

Next we derive a nonlinear equation for the global 
steady state of the laser, that includes a number of pulses 
and their accompanying pedestals, by applying the gain 
balance principle to the pulses, sidebands, and cw com- 
ponents of the waveform. We find that the sideband in- 
tensity and the pedestal width increase by a large fac- 
tor when the pumping is increased with a constant num- 
ber of pulses, and then decrease abruptly when another 
pulse is formed, so that the properties of single-pulse 
sidebands display an oscillatory, approximately periodic, 
dependence on the pump power. These theoretical re- 
sults are again favorably compared with experiments in 
a mode locked fiber laser. We conclude by considering 
the implication of our results on the nature of sidebands- 
mediated interactions. 

b. Theoretical model Our theoretical analysis is 
based on the mode-locking master equation model for 
soliton lasers [HI, [l9|, [2l| , where the dominant dynamical 
processes are chromatic dispersion and Kerr nonlincarity, 
whose coefficients can be nondimcnsionalizcd by an ap- 
propriate choice of units of time and power. The master 
equation then takes the form 

id z vb = -d\vb - 2|V>| V 

+ ip g {z)gW\ (^ + 7^) 

+ ip s (z)(s(\^\ 2 )-l)^ + T(z,t) (1) 
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The first two terms on the right-hand-side are the afore- 
mentioned dispersion and Kerr nonlinearity in soliton 
units. The next term models the saturable gain of the 
laser amplifier. g[ip\ is the overall gain coefficient, and 
the square brackets signify that g depends on the entire 
waveform ip(t) rather than the instantaneous value of tp 
only, and 7 is the coefficient of parabolic spectral filter- 
ing. We will assume that the the gain saturation is slow 
compared to the cavity round trip time so that g is 

tfl/2 

determined by the overall power P = J \ip\ 2 dt hi the 

-t R /2 

usual manner 



So 



(2) 



1 + P/Ps 

where go and P s are the small signal gain and the satu- 
ration power, respectively. 

The next term on the right-hand-side of the master 
equation models the fast saturable absorber with trans- 
missivity s(|?/>|) — I, where I is the small signal loss, and 
s(0) = by definition. We do not assume a particu- 
lar form for the transmissivity function other than that 
s(0) < I and that it increases linearly at zero, s'(0) > 0. 
We assume for simplicity that the saturable absorber ac- 
counts for all the losses in the cavity. The final term is a 
Gaussian white noise source with covariance 



(T*(z,t)r(z',t')) = 2T8(z-z')6(t-t'), 



(3) 



where the constant T is the rate of internal and injected 
noise power. As mentioned above, the noise is a signifi- 
cant factor in the determination of the steady state; we 
conjecture that it is also an essential ingredient in the 
inter-pulse interaction. The sidebands are formed by the 
spatial inhomogeneity of the gain and loss processes in 
the laser, described the by p g (z) and p s (z) (respectively) 
that are periodic functions of the cavity roundtrip length 
L, normalized to j- J Q L dzp g%s dz = 1. 

The dominance of the dispersive effects means that the 
gain and loss terms in Eq. ([1} are proportional to a small 
parameter, and that the noise term is proportional to an 
independent small parameter. In spite of their smallncss, 
the gain, loss and noise are the crucial terms for the mode 
locking phenomena discussed here. At the same time, 
they also perturb the pulse properties that are dominated 
by the dispersive terms; these small perturbation will be 
neglected. 

c. Single-pulse side bands We begin our analysis as- 
suming conditions under which there is a single pulse 
in the cavity with fixed parameters in the steady state. 
Since the dominant terms in the master equations are 
the dispersion and Kerr effect, the pulse waveform is ap- 
proximately that of a nonlinear Schrodinger (NLS) soli- 
ton. Solitons are defined by four parameters, amplitude 
a, frequency, timing and phase. The gain and loss terms 
in the master equation fix the frequency to zero, and we 
can set the timing and phase to zero by an appropriate 
choice of origin, so that the soliton waveform is 



The soliton waveform ip s is not an exact solution of 
the master equation because of the gain, loss and noise 
terms. We therefore look for a solution of the form 



ip(t, z) = ip p (t, z) + tp b (t, z) + ip c (t, z) 



(5) 



that in addition to the pulse waveform, consists of the 
sidebands waveform i\) b geneated by the cavity inhomo- 
geneities, and the continuum wave ip c generated by the 
cavity noise. The three waveform components have dif- 
ferent characteristic time scales: the sidebands are nar- 
row resonances whose temporal width is, as shown below, 
inversely proportional to the gain-loss small parameter. 
This width is large compared with the pulse width, but 
small compared with the cavity roundtrip time — the scale 
of the continuum. 

Both tpb and ip c have low peak power, and can therefore 
be analyzed by the linearized master equation, although 
their total power J dt\tfj b \ 2 , J dt\ip b \ 2 can be of the order 
P. The effect of noise is negligible on ip b , that therefore 
satisfies the equation 



+ Pg (z)g[ 1 p] (l + jd 2 ) (ip p + Tp b ) 
+ p p {z)(s{\^ p \ 2 )-l)(^ p + ^ b ) 



(6) 



The right-hand-side of Eq. ([6]) retains terms that are of 
higher order of smallness than ij) b ; these terms are in 
effect not negligible for \t\ 3> a where ip p is itself small, 
and play a crucial role in the shaping of the side bands, 
as shown below. 

The discrete modes of the real- linear Eq . (PdI) express 
small variations of the pulse parameters [20 . \22\ , while 
iphi that consists of radiation emitted by the pulse, is a 
linear combination 



V'b(M) 



2^ 



(a u (z)u u (t) + a u (z)*v u (t)*) (7) 



of the first component of scattering states of the linear 
operator L, 



I ={-{i + 91 )uj 2 -ia 2 )) lU " 



\={{ l -g 1 )u 2 + l a 2 )[ V J 



that acts on two-component wave functions as 

A B 



with 



B* A" 



(8) 
(9) 

(10) 



vj)p(t,z) = ascch(at)e l 



(4) 



A = (i + g^)d 2 - ia 2 + tea seek (at) + s(trsech (at)) 

and B = 2iasech 2 (at) . It will be argued below that we 
can use approximate the space-dependent linear operator 
acting on ip b in Eq. ^ by its space average L. Within this 
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approximation the coefficients defined in Eq. ([7]) evolve 
according to 

dzCtu, =(-(* + 9l)u 2 -ia 2 )a u + (ga^ + /3 gu )p g (z) 

+ (-la u + Psu)Ps(z) (11) 

where (3 g and (3 S are the expansion coefficients of the forc- 
ing terms g(l + "fd 2 )(ascch(at)) and (s(a 2 sech 2 (at)) — 
l)asech(at) (respectively). As usual, the expansion coef- 
ficients are extracted by inner product with the adjoint 
eigenfunctions defined by 



L^l»)={{i- gi )J + ia*){J + ai)[ U ») 



'. uz 



(-(i + 91 )uj 2 - ia 2 )) I 



(13) 



normalized so that 

dt(u u (t)*u u ,(t) + v^tyv^it)) = 2nS(uj-u>') (14) 
dt(v u (t)u^(t) +u u {t)v LJ ,(t)) = (15) 



Since by assumption the dispersive terms are dynam- 
ically dominant, the eigenfunctions of L are close to the 
eigenfunctions of the linearized NLS equation [23[ 



u u (t) = e 1 



2iuje 



(uj + ia) 2 



asech(at) + 



- 2 sech 2 (ai) 



(u + ia) 



v u (t) = e l 



a sech (at) 
(uj + ia) 2 



i \ (w + ia) 2 . . 
U u (t) = a- r^U-ui-t) 

(uj — %ay 

, . (uj + ia) 2 , 
v u (t) = -a- — jV- u (-t) 



(uj — iaf- 



(16) 
(17) 

(18) 

(19) 



so that 



= 7^ ( " + TZ I + ia 2 d 2 )asech9 



2n (uj — ia)' 



2iue at asech(at) + a sech (at) 



(uj + ia) 2 

Pslu — —-, — T^Xo. [ dt (s (a 2 sech 2 (at)) — l)ascch(at) 



2n (uj — ia)' 



(21) 



2iue at ascch(at) + a sech (at) 
(uj + ia) 2 

The solution of Eq. (JTTJl is 

a u (z) = J d z > e (-" 2 (i+9"/)-™ 2 )(z-*') ( 22 ) 

x e /;, dz"{gp g {z")-l P3 {*"))(p guPg ( z ) + p su] p s ( z )) 



As observed by Gordon and Kelly [TJ-Q , the amplitudes 
a u are driven resonantly if the nonlinear frequency shift 
a 2 + uj 2 is an integer multiple of the cavity based wave- 
length 2E; spectrally, therefore, the dispersive waves con- 
tain a discrete set of sidebands at frequencies ±w n , 



2- 



n = 1,2, 



(23) 



The nth sideband is forced mainly by the nth harmonic 
of the gain and loss modulation functions p g j, but also 
by the z dependence of the of the gain and loss terms in 
the exponent in Eq. (|22|) . The latter modulation is small, 
however, if we make the simplifying assumption that the 
total net loss per roundtrip is small, so that the z' in- 
tegration in Eq. (|22[) is effective over many roundtrips. 
This is the the usual assumption underlying the Haus 
master equation [l8| , and is also consistent with the rel- 
ative weakness gain and loss processes in the dynamics 
dominated by the chromatic dispersion and Kerr nonlin- 
carity that is studied here. 

In this case the mean net loss in the z" integration in 
Eq. (|2"2")l dominates over the variable parts of the gain and 
loss, and the latter may be neglected. The same assump- 
tion justifies the approximation of the time-dependent 
linear operator in Eq. (HU), the equation of motion for ipb, 
by the fixed operator L that was used to derive Eq. (jlll) . 

It now follows that the dispersive wave amplitudes are 

a u (z) = f dz ' e i-^ 2 +^)+ 9 {i-i^)-i){^') 

X (P g u,Pg(z) + P su p s (z)) (24) 

and defining the nth Fourier components p g ^ n and p S; „ 
of p g and p s (respectively), the coefficients of the nth 
sideband in the steady state are 



(Pg,nPgu + Ps,nPsu)e g 

I + g(~fuj 2 - 1) + i(uj 2 - uj 2 ) 



(25) 



We wish to characterize the sidebands by their ordi- 
nary Fourier spectrum, obtainable from Eq. (|7|). 



4>bn(u,z) = e m z J — (aw(z)w w ' w + a n ^(z)*v^ luj ) 

(26) 

where u u > u = J dte tuJt u u >(t), and is defined sim- 
ilarly. As a transform of a rapidly varying function, u 
is wideband — its uj bandwidth for fixed uj' is compara- 
ble with the soliton bandwidth — but it is also singular 
for uj = uj' . The smooth part of u generates a small 
deformation of the soliton waveform in Eq. (|2l)]) that is 
unimportant for the present purpose of characterizing the 
sideband spectrum. We therefore focus on the singular 
part iiuj'uj of Uu'u, that is determined by the ±oo asymp- 
totes of u u i (sec Eq. (jTHJ) ) to 



i^uj'uj — 2tt r i _ \o 



(uj + ia) : 



S(uj - uj') + 



4uja 



(uj + ia) 2 uj — uj 



7 (27) 
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where P denotes principal part, v^'u, = since t) u (t) 
tends to zero exponentially for t — > ±00. 

Now we can carry out the u/ integration in Eq. (|26[) 
and obtain the spectrum of the nth sideband 



I + g(jLJ 2 - 1) + i{uj 2 - ujD 



(28) 



The spectral width (half width at half maximum) of 
the sideband is given by 8uj = (I + g(juj 2 — l))/2w„ <C 1; 
temporally, therefore, the dispersive wave is a pedestal 
of frequency uj n with an exponentially decaying envelope 
centered at the pulse, whose decay time scale is much 
wider than the soliton width. 

Experimentally, the sidebands appear as a series of 
sharp peaks in the pulse spectrum on the background 
of the wide soliton spectrum, (see Fig. [TJ). A distinctive 
feature of the sideband spectrum (f2"5]l is the dependence 
of the overall phase of the sideband on the measurement 
position in the cavity; as a consequence, the sideband 
spectrum interferes constructively or destructively with 
the soliton spectrum, depending on the placement of the 
output coupler, exhibiting a sharp notch feature for some 
placements as the result of destructive interference. Figs. 
HH2] show the comparison between the measurement and 
calculation of the soliton and the first two sidebands spec- 
trum for two cavity positions. The main qualitative dis- 
crepancy between theory and experiment is the asymme- 
try between the left and right sidebands that is lacking 
in the theory. The most likely source of this discrepancy 
is the assumption of a nonchirped zero frequency pulse, 
that only approximately holds in the experiments. 
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FIG. 1: Measured spectra of the soliton and its sidebands, 
showing the notch formation. The two figures correspond 
to two measurement locations in the same cavity that are 
separated by one quarter of cavity length. These mesurements 
agree well with the theory. 

d. The continuum component The continuum ip c is 
the dominant component of the waveform for most of its 
temporal extent, where both i/j p (t) an d i>b(t) are negligi- 
ble. Therefore, the nonlinearity and interaction with the 
pulses are unimportant for its dynamics, and it is natu- 
ral to express it in terms of the ordinary Fourier modes 
i\) c (to) that satisfy 

d z ipc = -iu> 2 ij) c +pg(z)g (l - ijj c -p a (z)l'4> c +f(oj,z) 

(29) 

Because the noise and continuum extend throughout 
the entire roundtrip time, we must use discrete fre- 
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FIG. 2: Theoretical calculation of the pulse and sideband 
spectrum with: a = l,g = 0.1,7 = 0.01,2 = 0.11, s = 
0.1(H 2 - \ip 4 \),t R T = hn 3 ,L = 2.56,p 9 = e i7l/2 ,vs = e** . 
The cavity position in the left figure is z = L/5, and in the 
right is z — 0. 



quencies to label their Fourier transform. It then fol- 
lows from Eq. © that f(w, z) == J dte~ iult T{t, z) has 
zero mean and correlation function {T(oj, z)T*(co', z')) — 
2Tt R 8(z - z')6 u ,u>- The solution of Eq. is similar to 
(f2"2")l . and by the same arguments that lead to Eq. (|2"4"|) 
we again approximate the integrand in the exponent by 
its mean value. The resulting expression implies that 
(V'c(w)) = and 



$cM&(w')> 



l + g(juj 2 - 1) 



(30) 



e. Gain balance and multipulse sidebands In the 
analysis presented so far, the pulse parameters and the 
overall saturated gain g were assumed fixed and given. 
In the steady state these are variables, determined along 
with the number of pulses as a solution of the optical 
equation of motion Eq. (UJ). Since the equations of mo- 
tion are random, the result is a statistical steady state. In 
previous works (lol [l7j , statistical light-mode dynamics 
(SLD) theory was used to study this problem, and ap- 
plied to multipulse mode locking in 0, [H| . The global 
mode locking analysis determines in particular the prop- 
erties of the sidebands, and in this way allows us to reach 
our goal of describing the width and power of the side- 
bands as a function of the laser parameters. 

Here we study the statistical steady state by the gain 
balance method, deriving equations of motion for the 
power in the three components of the optical waveform, 
that is, the power P p in the pulse waveform that com- 
prises zero or more pulses, the sideband power and 
the continuum power P c . The three waveform compo- 
nents are characterized by well-separated time scales, and 
the total power can therefore be calculated as a sum of 
the powers of the individual components. One looks for 
a steady state where the three components are subject 
to the same saturated gain, and finally the gain itself is 
determined self-consistently from Eq. @. 

We will assume that the multipulse waveform is of 
the simplest kind [l2| that is also the most commonly 
observed experimentally, consisting of k pulses of equal 
amplitude a. In the equation of motion for the pulse 
amplitude we may neglect the noise term in the mas- 
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FIG. 3: The theoretical calculated values of the sidebands 
peak intensity, width, and total power Pi, along with the sat- 
urated gain coefficient g as a function of the small signal gain 
go- The same parameters as in Fig. [2] were used in the calcu- 
lation, except that 7 = 10~ 3 and Ps = 1. Two noise injection 
rate values are shown: T = 10 -3 in continuous black line, 
and T = 3 x 1CF 3 in dashed red line. The number of pulses 
is indicated. All values are given in natural soliton units. 



ter equation, and use the results of soliton perturbation 
theory @, to write 



da 



— =ReJ dttl>;(gp g (z){l + 7 fl d?Wp 

+ Ps (z)( s (IVvl 2 )-0Vv) 



= 2gp g (z)a(l - ja z )+p s (z){2s a (a) - 2l)a (31) 



where 



So (a) = — s(a 2 sech 2 (at))a 2 sech 2 (at)dt (32) 
2a J 

The pulse amplitude changes periodically during its prop- 
agation in the cavity; for our purposes we need the mean 
pulse amplitude a(z) = j- J a{z + z')dz' . Under the as- 
sumption of weak gain and loss processes the dynamics 
of a is obtained by the space average of Eq. (|31[). that 
gives in the steady state 



g (l-la 2 )-l + s a (a) = 



(33) 



This equation, together with P p = 2ka, determines the 
pulse part of the gain balance. 

Next, we calculate the power carried by the pedestals 
of the pulses. Each pulse generates a series of sidebands 
of the form given by Eq. (f2"5|) . The sidebands power is 
dominated by the leading, n = ±1 sidebands; since the 
different pulses that act as sources for the sidebands are 
not phase locked, we treat them as incoherent sources, 
and accordingly calculate the total sidebands power, as 





FIG. 4: The theoretical calculated values of the total pulse 
power P p and noise power P n as a function of the small signal 
gain go. The same parameters as in Fig. [3] The number of 
pulses is indicated. All values are given in natural soliton 
units. 



the sum of the individual sideband powers. The resulting 
expression is 



Ph = k- 



\Pg,lPgui +Ps,l/W| : 



2wi(J + fl (7w?-l)) 



(34) 



Finally the mean power of the the continuum is not 
changed by the presence of the sidebands, and is therefore 
given by the expression derived in (25j | 



Pr = 



Tt R 



VS7(* - .9) 



(35) 



smaller than the pulse power 2a, so that P c is of the same 
order of magnitude as P s . 

We now substitute P = P s + P b + P c into Eq. © and 
obtain a nonlinear equation for the saturated gain g that 
is easy to solve numerically. Once we know the value 
of g we obtain the steady state solution of the entire 
waveform. The steady state equation can have several 
solutions with different numbers of pulses k. In such cases 
the laser waveform can exist in several states [l2l [l3j , in 
a similar manner to the existence of metastable phases 
in thermodynamic systems undergoing first order phase 
transitions. The laser then exhibits hysteresis — its actual 
state depends on the history. 

Figs. [3H3] show the results of the calculation as the 
small signal gain go is varied. For very low gain the 
waveform is pure continuum, and its energy increases 
with the gain. When go is increased beyond a certain 
threshold, a pulse forms and along with it also a dis- 
persive wave. Because of gain saturation, the net gain 
must abruptly decrease, and along with it also the con- 
tinuum component. Further increase of go will mainly 
increase the sidebands and the continuum components, 
and slightly change the pulse amplitude, until the second 
threshold is met, and then the continuum as well as the 
sidebands power abruptly decrease again. This process 
then continues periodically when further pulse creation 
thresholds are reached. In addition, as the pumping is 
increased between pulse creation thresholds, the satu- 
rated gain increases, so that the net loss decreases and 
the sidebands become spectrally narrower, and accord- 
ingly temporally wider. 
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FIG. 5: The measured sideband spectrum characteristics: 
peak power, spectral width, and total energy, and the distance 
between adjacent pulses in states with two or more pulses. 

The theoretical predictions of the global sidebands 
characteristics agree well with experimental observation 
summarized in Fig. [51 made on a fiber ring laser mode 
locked by polarization rotation as described in [l3|. In 
both graphs the sideband width displays roughly peri- 
odic sawtooth behavior and the total sidebands energy a 
sawtooth behavior with an increasing linear bias. There 
is somewhat worse agreement in the peak intensity, but 
this quantity is sensitive to interferometric enhancement 
and reduction as explained above. 



In both theory and experiments the spectral width of 
the sidebands narrows with increasing pumping between 
pulse creation thresholds, implying a temporal widening 
of the pedestals. This widening shows a striking corre- 
lation with the increase of the inter-pulse distance, sug- 
gesting that the sidebands play a role in the inter-pulse 
interactions Hol|. 

/. Conclusions and outlook We presented a funda- 
mental theory of the formation of sidebands and their 
effect on the statistical steady state in multipulse mode- 
locked soliton lasers, that explains the main experimental 
observations, including the dependence of the sideband 
spectrum on the measurement position in the cavity, the 
growth of the sideband energy and temporal width when 
pumping is increased, and the abrupt attenuation of the 
sidebands when a new pulse is created in the cavity. 
We also found strong correlations between the temporal 
width of the sidebands and the spacing between adjacent 
pulses in pulse bunches, giving further evidence for the 
role of the sidebands in inter-pulse interaction. However, 
the most natural conclusion from our observation is that 
the sidebands generate repulsive interactions, and that 
additional attractive interactions are needed to explain 
the ubiquitous formation of pulse bunches. Moreover, 
the phase and timing jitter of the pulses leads us to con- 
jecture that the mechanism of interaction is incoherent. 
We postpone the in-depth study of these question to a 
future publication. 
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